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Abstract 

The critical behavior in the short-time dynamics for the random-bond Potts fer- 
romagnet in two-dimensions is investigated by short-time dynamic Monte Carlo simu- 
lations. The numerical calculations show that this dynamic approach can be applied 
efficiently to study the scaling characteristic, which is used to estimate the critical ex- 
ponents 0, f3/v and z, for the quenched disorered systems from the power-law behavior 
of the feth moments of magnetizations. 
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1 Introduction 



The understanding of the role played by quenched impurities on the nature of phase tran- 
sitions is one of the significant subjects in statistical physics, and it has been a topic of 
substantial interest for many authors in the last two decades |, ||, |, ||, |, 0, |, ||. Ac- 
cording to the Harris criterion fill], quenched randomness is a relevant perturbation at the 
second-order critical point when its specific-heat exponent a of the pure system is positive. 
Following the earlier work of Imry and Wortis[jl)| who argued that a quenched disorder could 
produce rounding of a /irst-order phase transition and thus induce second- order ^\y&s,^ tran- 
sitions, the introduction of randomness to systems undergoing a first-order phase transition 
is comprehensively considered. It was shown by Hui and Berker that the bond randomness 
can have a drastic effect on the nature of a first-order phase transition by phenomenological 
renormalization group arguments 0, and the feature has been placed on a firmer basis with 
the rigorous proof of vanishing of the latent heat 0. Their theory was numerically checked 
with the Monte Carlo (MC) method by Chen, Ferrrngberg and Landau (CFL) ^, where 
the eight-state Potts model was studied with random-bond disorders. Experimental evi- 
dence in two-dimensional systems was found that in the order- disorder phase transitions of 
absorbed atomic layers, the critical exponents are modified, in the addition of disorder, from 
the original four-state Potts model universality class in the pure case [l^ . On the other 
hand, no modification is found when the pure system belongs to the Ising Universality class 



T^. The theoretical study of such kinds of disordered systems is also an active field where 



a resort to intensive MC simulations is often helpful [|T^, |T^, |T8|, |T9[ . 

It is well known that the pure Potts model in two-dimensions (2D) has a second order 
phase transition when the number of Potts state g < 4 and is first order for g > 4. As the 
specific heat exponent a of the pure system is always positive for g > 2, the disorders will be 
the relevant perturbation for the Potts model. As a result, all the transitions are second order 
for all 2D g-state Potts models in the presence of quenched disorders, and the impurities 
have particularly strong effect for g > 4, even changing the order of the transitions. 

In this paper, we discuss the dynamic scaling features of random-bond Potts model 
(RBPM) through MC simulations, to estimate the critical exponents. We consider the 
important questions of whether there exists an Ising-like universality class for the RBPM 
and how is the critical behavior affected by the introduction of disorder into the pure system 



T7| . The large-scale MC results by the CFL and in ref. [0 suggest that, in 2D, any random 
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systems should belong to the pure Ising universality class. These results are also coherent 
with experiment ||12[. In recent papers]^, ||, however, Cardy and Jacobsen studied the 
random-bond Potts models for several values of q with a different approach based on the 
connectivity transfer matrix (TM) formalism of Blote and Nightingale |^D|] Their estimates 
of the critical exponents led to a continuous variation oi (3/v with g, which is in sharp 
disaggreement with the MC results for g = 8 |@, ||. We hope that the resulting critical 
behavior measured in this paper will play a role in settling this controversy. Furthermore we 
will test the short-time dynamic (STD) MC approach for the first time, to study the spin 
systems with the quenched disorder and to show its efficiency by numerical studies, which is 
also one of our main aims of this paper. 



2 Model and Method 

The Hamiltonian of g-state Potts model with quenched random interactions can be written, 

-(3H=Y. > >o, (1) 

<i,i> 

where the spin a can take the values 1,- ■ ■ q, (3 = l/ksT the inverse temperature, 6 is the 
Kronecker delta function, and the sum is over all nearest-neighbor pairs on a 2D lattice. The 
dimensionless couplings Kij are selected from two positive (ferromagnetic) values of Ki and 
K2 = rKi, with a strong to weak coupling ratio r = K2/K1 called as disorder amplitude, 
according to a bimodal distribution, 

P{K)=pS{K-Ki) + {l-p)S{K~K2) . (2) 

When p = 0.5, the system is self-dual and the exact critical point can be determined by | |21|| , 

(e^^-l)(e^— l) = g. (3) 

where Kc and K'^ are the corresponding critical values of Ki and K2 respectively at the 
transition point. While r = 1 corresponds to the pure case and the critical point is located 
at Kc = log(l + ^/q~) and the phase transitions are first-order for g > 4. With additional 
random-bond distribution, however, new second-order phase transitions are induced for any 
of g-state Potts models and the new critical points are determined according to Eq.(^) for 
different values of disorder amplitude r and state parameter g. 



2 



In this work we chose q = 8, which is known to have a strong first-order phase transition, 
in the hope that we would find a new second-order phase transition caused by the quenched 
disorder to show the effect of impurities on the first-order systems. The strength of the 
disorder was chosen for several values of r, as was done in ^, |TD[, to check the Ising-like 



universality class. To minimize the number of bond configurations needed for the disorder 
averaging, we confined our study to the bond distributions in which there are the same 
number of strong and weak bonds in each of the two lattice directions. This procedure 
should reduce the variation between different bond configurations, with no loss of generality. 

We performed our simulations by the STD method on the 2D square lattices with pe- 
riodic boundary conditions. This dynamic MC simulations have been successfully performed 
to estimate the critical temperatures Tc and the critical exponents 6, (3, v and dynamic expo- 
nent z for the 2D Ising model p3| and the 2D 3-state Potts model , since for both models 



there exist second order phase transitions. Recently this approach has also been extensively 
applied for the Fully Frustruated XY model and spin glass systems to study the critical 
scaling characteristics to estimate all the dynamic and static critical exponents [^, |2^ . 
Traditionally it was believed that universality and scaling relations can be found only in 



the equilibrium stage or long-time regime. In Ref.|]5D|, however it was discovered that for a 
magnetic system in states with a very high temperature T ^ which is suddenly quenched 
to the critical temperature and then evolve according to a dynamics of model A |^], there 



emerges a universal dynamic scaling behavior already within the short-time regime, which 
satisfies, 

M^'^) {t,T,L,mo)= h-^^'^M^^^ {b-H, b'/'r, 6" ^L, 6"°mo) , (4) 

where M^''^ is the fcth moment of the magnetization, r = (T — Tc)/Tc is the reduced tem- 
perature, (3 and v are the well known static critical exponents and 6 is a scaling factor. The 
variable xq, a new independent ex^'^newi, is the scaling dimension of initial magnetization mo. 
This dynamic scaling form is generalized from finite size scaling in the equilibrium stages . 
Importantly the scaling behavior of Eq.(^ can be applied to both dynamic exponent mea- 
surements and the estimates of the static exponents originally defined in equilibrium. 

We begin our study on the evolutions of magnetization in the initial stage of the dynamic 
relaxation starting at very high temperature and small magnetization (mo ~ 0). For a 
sufficiently large lattice (L oo), from Eq.(§) by setting r = 0, 6 = t^/^, it is easy to derive 
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that 

When k = 1 we get the most important scahng relation on which our measurements of the 
critical exponent 9 are based, 

M{t) ~ mo^^ = {xo- P/iy)/z. (6) 

As a result, the magnetization undergoes an initial increase at the critical point Kc after a 
microscopic time tmic- This prediction is supported by a number of MC investigations which 



have been applied to detect all the static and dynamic critical exponents p3| , as well as 



the critical temperatures p5| , p8| . The advantage of the dynamic MC simulations is that it 
may eliminate critical slowing down since the measurements are performed in the early time 
stages of the evolution where the spatial and time correlation lengths are small. 

In our simulations, the time evolution of M{t) is calculated through the definition 

Here Mq = max(Mi, M2, ■ ■ ■ , Mg) with Mj being the number of spins in the ith state among 
q states. < ■ ■ ■ > denotes the initial configuration averages over independent random num- 
ber sequences, and [■ ■ ■] the disorder configuration averages over quenched random-bond 
distributions. = is a number of spins on this square lattice and g = 8 is chosen. 

The susceptibility plays an important role in the equibrium. Its finite size behavior is 
often used to determine the critical temperature and the critical exponents 7/z/ and P/u 
For the STD approach, the time-dependent susceptibility (the second moment of the 
magnetization) is also interesting and important. For the random-bond Potts model, the 
second moment of the magnetization is usually defined as 

M(2)(t) = ![(< M\t) >-< M{t) >2)] . (8) 

To study the scaling behavior of the second moment of magnetization, we have to take 
the initial states of mo = to start the relaxation processes. Because the spatial correlation 
length in the beginning of the relaxation is small compared with the lattice size L'^ in the 
short-time regime of the dynamic evolition, the second moment behaves as M^^^(t, L) ~ L~'^. 
Then the finite size scaling Eq.(^ induces a power-law behavior at the critical temperature, 

M(2) {t)r~.t\ y = {d- 2P/u)/z. (9) 
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From the scaling analysis of the spatial correlation function we easily realize the non- 
equibrium spatial correlation length ^ ~ t^^^. Therefore M^'^\t) ~ ^('^-2/3/i^)_ 

In the above considerations the dynamic relaxation process was assumed to start from a 
disordered state or with small magnetization mg. Another interesting and important process 
is the dynamic relaxation from a completely ordered state. The initial magnatization locate 
exactly at its fixed point mo = 1, where scaling of the form, 

M^^\t, r, L) = b-'"^/''M^''\b-H, b^/^T, b-^L), (10) 

is expected. This scaling form looks to be the same as the dynamic scaling one in the 
long-time regime, however, it is now assumed already valid in the macroscopic short-time 
regime. 

For the magnetization itself, b = t^^^ yields, for a sufficiently large lattice, 

M(t, r) = t-^/^W(l, t-^/^V). (11) 

This leads to a power-law decay behavior of 

M{t,T)=r^\ ci = p/vz, (12) 

at the critical point (r = 0). The formula can be used to calculate the critical exponents (3/v 
and z. For a small but nonzero r, the power-law behavior will be modified by the scaling 
function M(l, t~'^/'^^r), which has been used to determine the critical temperatures p8| , ^ . 
Furthermore, by introducing a Binder cumulant 

Uit L) - - 1 ri3) 

(M(t,L))2 ^' 



a similar power-law behavior at the critical point induced from the scaling Eq.([10D shows 
that, 

U{t, L) ~ t'^ , C2 = d/z , (14) 

on a large enough lattice. Here, unlike the relaxation from the disordered state, the fiuctu- 
ations caused by the initial configurations are much smaller. In pratical simulations, these 
measurements of the critical exponents and critical temperature are better in quality than 
those from the realization process starting from disordered states. 
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Table 1: The tendency and measured values of 6' as a function of the disorder amplitude r 
for different initial mo at the critical points Kc on the lattice 64^ lattice. 



mo 0.06 0.04 0.02 0.01 


e 




r = 2 0.310(8) 0.338(8) 0.350(8) 0.352(7) 


0.353(6) 


0.920185271... 


r = 5 0.160(6) 0.215(6) 0.252(5) 0.257(4) 


0.262(4) 


0.512307010... 


r = 8 0.106(5) 0.167(4) 0.208(4) 0.218(3) 


0.221(3) 


0.367963156... 


r = 10 0.090(4) 0.146(3) 0.193(3) 0.202(3) 


0.203(3) 


0.312655667... 



3 MC Simulations and Results 

As it was pointed out that the Heat-bath algorithm is more efficient than the Metropolis 



algorithm in the STD |2^, and the universality is satisfied for different algorithms, we only 
perform the MC simulations with the Heat-bath algorithm at the critical points of 2D eight- 
state RBPM for an optimal disorder amplitude r* = 10 which is located on the random fixed 
point regime with a largest value of central charge c = 1.5300(5) |jlO[- Samples for averages 
are taken both for over 300 disorder distribution configurations and about ~500 independent 
initial configurations on the square lattices with L up to 128. Statistical errors are simply 
estimated by performing three groups of averages with different random seed selectes for the 
initial configurations. It should be noted that, except for M(t), the measurements of M^'^\t) 
and U{t) are restricted to the initial states with itiq = or mo = 1. Importantly, it was 
verified that the critical exponents save the same value the same as those in the equilibrium 



or long-time stage of the relaxation [^. Therefore we can measure these exponents based 
on the corresponding scaling relation in the initial stages of the relaxation. 

We start our simulations to verify the power-law scaling behavior of M{t) with several 
values of disorder amplitude at the critical points Kc{r) (as shown in Table 1). The initial 
configurations are prepared with small magnetizations mo = 0.06, 0.04, 0.02 and exact zero 
states. In Fig.l, the time evolutions of the magnetization M{t) versus the disorder amplitude 
r on a 64^ lattice are displayed with a double-log scale. We can easily find that all the curves 
exhibit the power-law behavior predicted by Eq.(^). Thus 6 can be estimated from the 
slopes of the curves. The values of 6' as a function of the disorder amplitude r for small 
initial magnetization mo are presented in the Table 1. 

We then set mo = to measure the second moment of magnetization. The power-law 
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Table 2: The values of scaling exponents for the 2D q=8 RBPM with r = 10, measured 
from the scaling functions of M(t), M^'^\t) and U{t) respectively starting from both the 
random initial states and ordered states. Also listed are those for the 2D Ising and g = 3 
Potts models, and the 3D Ising model [E^, pi, . 



exponent 


mo 


2D RBPM 2D Ising 2D Potts 


3D Ising 


e 

y=id-2/3/u)/z 


~ 0.0 


0.197(4) 0.191(1) 0.075(3) 
0.438(6) 0.817(7) 0.788(1) 


0.108(2) 


ci = P/uz 
C2 = d/z 


= 1.0 


0.0390(6) 0.056(1) 0.065(1) 
0.518(9) 0.926(8) 0.934(9) 


0.2533(7) 
0.1462(12) 


2(3 /u = d-yz 
2(3 /v (exact) 




0.302(6) 0.240(36) 0.269(7) 
1/4 4/15 


1.034(4) 



behaviour of the second moment M'-^^(t) is observed in Fig. 2, where the curves for different 
lattice sizes are plotted. Again, they present a very nice power-law increase. Values of scaling 
dimension y = {d — 2(3 /v) / z determined from slopes of the curves during t = [10,200] are 
listed in the Table 2. 

We furthermore set mo = 1 to observe the evolution of the magnetization and the Binder 
cumulant, both should show the power-law behavior as predicted by the Eg. ([T^) and Eg. (p!3|) . 
Their curves are plotted in the Figs. 3 and 4 respectively. The values of the scaling dimension 
Ci = (3/vz and C2 = dj z are then estimated, as are also presented in Table 2. Now the results 
oiy^c\ = (31 vz and C2 = dj z can be used to estimate critical exponent (3lv shown in Table 2. 
For comparison, also listed in Table 2 are the correspoinding results of the scaling dimension 
for the Ising and g = 3 Potts models on 2D (3D) sguare (cubic) lattices, and in Table 3 we 
summarize the results of the critical exponent /9/z/ up to the present. 

4 Summary and Conclusion 

In this paper we have investigated the short-time critical dynamics of the random-bond Potts 
model on 2D lattices to verify whether it has a second order phase transition in the Ising-like 
universality class, by a STD study. Dynamic scaling haviour was found, and has been used 
to estimate the critical exponents 6', z and (3lv. Our main results are summarized in Table 
2, which are obtained from the slopes of power-law curves of M(t), M(2)(t) and U{t) in the 
double-log scales by least-sguare fits. 
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Table 3: The results of magnetic scaling exponent j3/v estimated by different methods for 
the 2D eight-state RBPM. 



Authors 


r 




Technique 


CFLi 


2 


0.118(2) 


MC 


Cardy and Jacobsen[|] 


2 


0.142(4) 


TM 


Chatelain and BercheQ 


10 


0.153(3) 


MC 


Piccog 


10 


0.153(1) 


MC 


Present work 


10 


0.151(3) 


STD 



Our work shows that for the RBPM, there exists a power-law behavior, which is the 
typical feature of a continuous phase transition in the STD processes. The r-dependence 
characteristics of the dynamic exponent 9 gives evidence that their dynamic MC behavior 
is different from the pure Ising model, and the values of magnetic exponent (3/v m our 
calculation seems to be the same as that given by the TM formula [Q, but not as that by 
the CFL IQ, . Furthermore we found that the values of dynamic exponent 9 depend on the 
disorder amplitude r, and it suggests that the dynamic exponent z may also depend on the 
r, which would be interesting to examine for future study. 

In conclusion, this study presents numerical evidence that the quenched impurities in 
the RBPM can induce new second-order phase transitions, but they appear not always 
to belong to Ising-like universality class, although the result of new critical exponent 9 is 
the same for both the r = 10 RBPM and the Ising model within errorbars by present 
calculations. Second, as the effect of critical slowing down in the equilibrium stage for the 
RBPM is more severe than that for the pure systems, the cluster algorithms, up to now, 
have been frequently used in MC simulations of the 2D RBPM ||, ^ |I0|, |T9[. In present 
paper, however, we have applied the STD MC simulations which use local updating schemes 
for the 2D RBPM, and the fact that dynamic MC simulations can avoid the critical slowing 
down in the STD processes, where the spatial correlation length is still small, makes it easier 
to calculate the critical exponents. An important subject for further study is, for example, 
how the dependence of dynamic critical exponent l3/v on the state parameter q and the 
disorder amplitude r by a systematic simulation using the STD method in order to clarify 



the crossover behaviour from the random fixed point to a percolation-like limit ||33[ . This is 
being studied at present. 
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FIGURE CAPTIONS 



FIGURE 1: The time evolution of magnetization showing the r-dependence of 6, plotted 
with a double-log scale on a lattice of 64 x 64 with = 0.01. 

FIGURE 2: The time evolution of second moment of magnetization starting from ab- 
solute random states, plotted with a double- log scale on lattices of 32 x 32, 64 x 64 and 
128 X 128. 

FIGURE 3: The power-law decay of magnetization starting from fully ordered states, 
plotted with a double-log scale on lattices of 32 x 32, 64 x 64 and 128 x 128. The finite size 
effect is obvious when the lattice size L < 64. 

FIGURE 4: The time evolution of Binder comulant starting from fully ordered states, 
plotted with a double-log scale on lattices of 32 x 32, 64 x 64 and 128 x 128. 
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Fig.1 : Evolution of magnetization in a log-log scale with m_0=0, 
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Fig. 2: Evolution of susceptibility from random-initial states 




10-^ L 
10 



100 

t (MC sweep) 



Fig. 3: Power-law decay of magnetization from ordered-initial states 
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Fig. 4: Evolution of Binder cumulant from ordered-initial states 
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